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Abstract. We report the results of a systematic compar- 
ison between the vertically averaged model and the ver- 
tically explicit model of steady state, Keplerian, optically 
thick a-discs. The simulations have concerned discs cur- 
rently found in three different systems: dwarf novae, young 
stellar objects and active galactic nuclei. In each case, we 
have explored four decades of accretion rates and almost 
the whole disc area (except the narrow region where the 
vertically averaged model has degenerate solutions). We 
find that the one layer approach gives a remarkably good 
estimate of the main physical quantities in the disc, and 
specially the temperature at the equatorial plane which is 
accurate to within 30 % for cases considered. The major 
deviations (by a factor < 4) are observed on the disc half- 
thickness. The sensitivity of the results to the a-parameter 
value has been tested for 0.001 < a < 0.1 and appears 
to be weak. This study suggests that the "precision" of 
the vertically averaged model which is easy to implement 
should be sufficient in practice for many astrophysical ap- 
plications. 
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1. Introduction 

Recent investigations on MHD turbulence have confirmed 
the existence of a mechanism able to extract angular mo- 
mentum from weakly magnetized accretion discs (Bal- 
bus & Hawley, 1991; Fleming, Stone & Hawley, 2000; 
Miller & Stone, 2000). Many aspects of accretion are, how- 
ever, still not understood. In the meanwhile, present disc 
models mostly incorporate the a-prescription (Shakura & 
Sunyaev, 1973) and are of increasing complexity. For in- 
stance, models of steady state, Keplerian, vertically strat- 
ified a-discs (e.g. Cannizzo, 1992; Hameury et al., 1998; 
DAlessio et al., 1999) should definitely supplant the more 
widespread, vertically averaged version (Pringle, 1981) 
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which, despite a great simplicity and flexibility, cannot 
account for fine physical effects, remains in essence some- 
what limited and contains some artefacts, like degener- 
ate solutions (CoUin-Souffrin & Dumont, 1990). Vertically 
explicit disc models offer in principle a much better reli- 
ability but, due to the uncertainty on the a-formalism, 
they require an extra prescription specifying the depth- 
dependence of the viscous energy release (Meyer & Meyer- 
Hofmeister, 1982; Mineshige & Osaki, 1983; Cannizzo & 
Cameron, 1988; Hameury et al., 1998). From a numeri- 
cal point of view, the bi-dimensional problem which in- 
volves partial/ordinary differential equations is much less 
trivial to solve than the vertically averaged problem, also 
meaning a much larger computational time (by a factor 
^ 100 — 1000 say, corresponding to the typical resolu- 
tion in the ^-direction). It is therefore legitimate to ask, 
from a quantitative point of view at least, which benefits 
the vertically explicit model really brings. The aim of this 
report is to show how, for steady Keplerian a-discs, the 
one layer approach and the vertically explicit approach 
compare. The main motivation is both to seek where the 
differences appear in an a-disc and to quantify these as 
functions of the accretion rate and value of the viscosity 
parameter. This investigation concerns three generic sys- 
tems containing a disc: an active galactic nucleus, a young 
stellar object and a dwarf nova. 

2. Vertically explicit vs. vertically averaged 
models: hypothesis and restrictions 

We compare the structure of steady state Keplerian accre- 
tion discs computed from the vertically averaged model 
(hereafter, the VA-model) and from the vertically explicit 
model (hereafter, the VE-model). The difference between 
the two resides therefore in the resolution of the vertical 
stratification. This analysis is carried out in the frame- 
work of the a-theory (Shakura & Sunyaev, 1973), with a 
uniform a-parameter value. For a detailed description of 
the models, the reader is referred to Hure (1998, 2000). 

In order to give a meaning to the comparison, the two 
models are made as close as possible and reduced to their 
simplest form. In particular, we use the same Rosseland 
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Fig. 1. Midplane temperature (top) and disc half-thickness (bottom) versus the radius for M = 0.61 M0, M ~ 10 
Mg/yr and a = 10"^ (DN; left), M = 1 Mq, M = lO"'^ Mg/yr and a = IQ-^ (YSO; middle) and Af = 10^ Mg, 
M = 10~^ Mg/yr and a = 0.1 (AGN; right). The systematic comparisons between the VA-model [dashed lines) and 
the VE-model {solid lines) concern the hot and cold branches, outside the zone where the VA-model has a multivalued 
solution (MVS). 



and Planck opacity grids and equation of state correspond- 
ing to a solar mix at LTE for both models. Convective 
transport is not taken into account in the VE-model, for 
several reasons. First, it cannot be included in the one 
layer version. Second, it plays a role in regions of high 
opacity, namely where hydrogen recombines. In the frame- 
work of the Mixing Length Theory (e.g. Cox & Giuli, 
1968), the disc region mainly involved has a small size and 
convection does not drastically change the disc structure 
(Meyer & Meyer- Hofmeister, 1982; Cannizzo, 1993). Note 
that the weak effect of convection is sustained in the work 
by Gu, Vishniac & Cannizzo (1999). Third, as explained 
below, we shall not analyze the region where the one layer 
has multiple solutions and which precisely encompasses 
the convectively unstable zone. 

Also left aside is external irradiation. It depends on 
many parameters and has little effect on the global quan- 
tities (midplane temperature, surface density, etc.) as long 
as the disc remains optically thick. Regarding the depth- 
dependent viscosity prescription, we adopt the most com- 
mon parameterization of the heat flux gradient, namely 



dz 2 



(1) 



where P is here the total pressure and J^k is the Keple- 
rian angular velocity. We also neglect the disc self-gravity. 



3. Results 

The disc structure depends mainly on three input param- 
eters: the central mass M, the accretion rate M and value 
of the a-viscosity parameter. The departures between the 
VA-model and the VE-model are measured at the same 
radius R by the ratio 

^VE — model 

^{x) — — TFI TT (2) 

^ ^ J.VA— model ^ ' 

where a; is one of the following quantities: the temper- 
ature Tmid (the subscript "mid" refers to the midplane), 
half the disc thickness h (taken as the altitude of the 
photosphere's base in the VE-model and as the pressure 
scale height in the VA-model), the density Pmid, the total 
surface density St = P'^^ ^^"^ the optical thickness 
Tt = 2 /q^°° Kpdz. Note that e is a function of the radius. 

Figure |l| displays h{R) and Tmid(^) obtained from the 
two models for parameter values (M, M, a) typical of a 
dwarf nova (DN) disc, a young stellar object (YSO) disc 
and a disc in an active galactic nucleus (AGN). We see 
that the VE-model always gives a slightly thicker and hot- 
ter disc than the VA-model, but the agreement is globally 
rather good. The other quantities, not shown here, com- 
pare very well too. 

In order to test the sensitivity of these results to the 
values of the input parameter values, we have performed 
a series of computations by varying the accretion rate and 
the value of a, for any of the above three systems. This 
amounts to about 29000 runs in total. Since the one layer 
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system 


central mass (in AIq) 


Q-parameter 


radii 


accretion rate (in Mq/jx) 


number of runs 


AGN 




0.01 and 0.1 


20 - lO^i^s 


10"* - 0.4 


9372 


DN 


0.61 


0.01 and 0.1 


10^ - 10" cm 


10-1^ - 10"** 


7742 


YSO 


1 


0.001 and 0.01 


0.05 - 50 UA 


lO"'' - lo-'^ 


11534 



Table 1. Domain of calculation, depending on the system (i?s is the Schwarzschild radius of the black hole). In all 
cases, the (log A/,logi?)-plane is explored with constant logarithmic steps corresponding to = ^ and = for 
the accretion rate and radius respectively. The last column gives the number of runs performed for the comparisons. 




-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 

LOG[Accrotion rata in M„yr'] 



Fig. 2. Iso- values of e(Tmid) in the (log i?, log A/)— plane, 
in the AGN case (M = 10* and a = 0.01). The MVS 
region (see Fig. I and Eq.(|)) is shown in grey. The dashed 
line labeled (3 = \ marks the limit between the gas pres- 
sure supported disc (at large radii) and the radiation pres- 
sure supported disc. 

model yields a multivalued solution (MVS), not present 
in the VE-model (Cannizzo, 1992), the comparison has 
been restricted to the hot, inner and cold, outer branches 
located on both sides (see Fig. |l|). The approximate loca- 
tion of the MVS region will be specified later. 

Figure || displays iso-values of the quantity e(T,„id) 
in the (log M, log i?)— plane in the AGN case, for the 
canonical values M = 10* Mq and a = 0.1. From this 
data grid, the maximum and minimum deviations on the 
temperature ratio can easily be extracted. So, we find 
1 < e(rinid) ^ 1-2 over the radial domain and for accretion 
rates indicated in Table |l|. We have repeated this proce- 
dure for the other four key quantities /i, Pmidj "^t and Tt, 
and for a = 0.01 as well. Table || summarizes the results. 
We notice that the agreement between the two approaches 
is excellent, indeed. The midplane temperature is slightly 
underestimated with the VA-model, by a factor never ex- 
ceeding 20% (in the cold outer disc) for any a, which is 
remarkable given the crudeness of this model. The other 
quantities are a little bit more affected by the vertical av- 
eraging: we have a mean factor w 2 on the central density 
and K, 1.5 on the optical thickness and surface density. 
The largest departures are observed on the disc thickness. 



a = 0.1 




hot, inner branch 


cold, outer branch 


e(r,nid) 


1.0 - 1.2 


1.1 - 1.2 


e{h) 


1.2 - 2.8 


2.2 - 3.9 


e(pmid) 


0.8-3 


1.0 - 2.8 


e(St) 


1.1 - 2.2 


1.3 - 1.6 


<n) 


1.0-2.3 


0.6 - 2.4 


a = 0.01 




hot, inner branch 


cold, outer branch 


e(rinid) 


1.0 - 1.2 


1.1 - 1.2 


e{h) 


1.8-3.8 


2.5-4.0 


e(Pmid) 


0.8-1.5 


1.0 - 1.3 




0.9-1.7 


1.0 - 1.3 




0.9-1.6 


0.6 - 1.6 



Table 2. Extreme values of e for an AGN disc, and for 
two different values of the a-parameter, outside the MVS 
region (see Fig. and Table ^ for the domain of compu- 
tation). 



Q = 0.01 




hot, inner branch 


cold, outer branch 


e(rmid) 


0.9-1.1 


1.0 - 1.3 


e{h) 


1.2 - 1.9 


1.0-3.3 


e(pmid) 


1.2 - 1.6 


0.4 - 1.4 


^(St) 


1.7-3.4 


0.9 - 1.8 




1.4 - 2.1 


0.3 - 1.9 


a = 0.001 




hot, inner branch 


cold, outer branch 


e(r,nid) 


0.9-1.1 


1.0 - 1.2 


e{h) 


1.1 - 2.6 


1.5-4.0 


e(Pmid) 


1.1 - 2.1 


1.0 - 1.4 


^(SO 


1.5 - 2.9 


1.3 - 1.8 




1.3 - 2.4 


0.7-2.0 



Table 3. Same as for Table ^ but in the YSO case. 



with a factor reaching « 4 in the cold, outer disc. The sen- 
sitivity to the value of the a- viscosity parameter is weak, 
even null for the temperature. However, the trend is that 
the deviations seem smaller with a smaller value of the 
a-parameter, except for h. 

We have carried out similar computations in the YSO 
and DN cases. The domains of calculation, specific to each 
system (Gullbring et al., 1998; Horne, 1998), are listed in 
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Q = 0.1 




hot, inner branch 


cold, outer branch 


e(Tmid) 


1.1 


1.1 - 1.7 


e{h) 


1.0-2.8 


1.8-3.8 


e(/5mid) 


1.1 - 2.6 


0.3-1.3 




1.5-3.2 


0.6-1.7 




1.3-2.9 


0.8-6.0 


a = 0.01 




hot, inner branch 


cold, outer branch 


e(Tmid) 


1.1 


1.1 


e{h) 


1.2 - 2.9 


2.6 - 2.9 


e(Pmid) 


1.2 - 1.4 


0.9-1.6 




1.5 - 1.8 


1.4 - 1.5 




1.3-1.8 


1.2 - 1.8 



Table 4. Same as for Table || but in DN case. 



Table |^ and the results are summarized in Tables || and ^ 
respectively. Globally, the conclusions drawn in the AGN 
case still hold. The best agreement between the two mod- 
els is observed in the DN case. In particular, the temper- 
ature ratio is constant (« 1.1) throughout, whatever the 
value of a. However, for a — 0.1 and for the lowest accre- 
tion rates examined, the optical thickness becomes close 
to unity in the outer disc and so the optically thick ap- 
proximation underlying the two models becomes doubtful. 
This produces a specially large ratio e{T^id) 1.7. 

We point out that the deviations between the two mod- 
els are largest in the vicinity of the MVS region, on the 
hot branch side (see Fig. |l|). Since our radial resolution 
(set to = 0.1 in the present calculations) is finite, it 
is expected that the "real" , maximum deviations (i.e. the 
deviations resulting from an infinite resolution) are prob- 
ably slightly different to those reported in the tables. 

Gathering the results obtained for the three systems 
considered here, we find that the MVS region which more 
or less encompasses the thermally unstable regime where 
hydrogen recombines (Meyer & Meyer-Hofmeister, 1982) 
is located at 

i?MVS - (1 ± 0.6) X 10i5a-°-2AfO-34MO-44 cm (3) 

where the mass Mq and accretion rate Mq are ex- 
pressed in Mq and MQ/yr respectively. The universal 
character of this expression can easily be understood from 
power law solutions of the VA-model (Hure, 1998), given 
that the MVS region (see Fig. ||) is roughly bounded by 
isotherms. 

4. Concluding remarks 

Given the ranges of parameter values (masses, accretion 
rates and viscosity parameter values) considered in this 
study, we can conclude that the vertically averaged model 
of steady state Keplerian a-discs is remarkably reliable, 
since it provides solutions very close to those obtained 



from the vertically explicit version. The midplane temper- 
ature appears as the quantity least sensitive to the vertical 
averaging (with a 30% deviation only) whereas the disc 
thickness seems the most affected (a factor ~ 4 in some 
cases). Problems dealing with the disc irradiation (even 
with a global energy budget) should therefore account for 
vertical stratification since the disc response is sensitive 
to the local flaring angle and radial gradient (Kenyon & 
Hartmann, 1987). 

This investigation demonstrates that the vertically av- 
eraged approach should be sufficient for many applications 
(e.g. Pringle, 1997; Burderi, King & Szuszkiewicz, 1998; 
Collin & Hure, 1999), in particular if the level of accuracy 
required on the disc structure is low or moderate. Note 
however that the tables given here can be used to "cor- 
rect" the vertically averaged solution if a better precision 
is necessary. More basically, these tables show the typical 
error on the knowledge of the disc structure, that is, some 
methodology effect. It is interesting, but comforting, to 
see that the model is much more sensitive to the physics 
it contains (like turbulent viscosity, convection, opacities, 
etc.) than to the method used to derive the disc structure. 

Acknowledgements. I thank M. Mouchet for reading the 
manuscript, D. Richard, and the referee, Dr John K. Cannizzo, 
for interesting remarks, specially regarding convection. 
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